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An alternative method is applied to the study of nucleon-nucleon(A'^A^) scattering phase shifts in 
the framework of extended quark delocalization, color- screening model(QDCSM), where the one- 
pion-exchange(OPE) with short-range cutoff is included. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is generally be- 
lieved to be the fundamental theory of the strong in- 
teraction. High energy processes are calculable due to 
the asymptotic freedom of QCD. Low energy processes 
are however difficult to be calculated directly from QCD, 
especially for multiquark systems, due to the infrared 
confinement and complexity. At present and even in the 
future, the QCD-inspired model will be a useful tool to 
explore strong interaction physics at low energy region. 

To study the baryon-baryon interaction and multi- 
quark system, the most popular model is the constituent 
quark model(CQM). There are different versions of CQM 
with different effective degrees of freedom. The Glashow- 
Isgur model (GIM) where the effective degrees of 

freedom are constituent quarks and gluons, describes the 
properties of hadrons successfully, and has also been ex- 
tended to the study of the NN interaction. The NN 
short-range repulsive core is successfully reproduced, but 
the intermediate range attraction is missing. To remedy 
this shortcoming, two approaches are developed. One is 
to invoke meson exchange again, this leads to the 'hy- 
brid' model d, 0, [1] , another approach is to extend the 
GIM, this leads to the quark delocalization color screen- 
ing modcl(QDCSM) 0,0]. Two approaches both fit the 
existed experimental data of the NN interaction, but 
only the QDCSM explains the long standing fact that the 
nuclear force and molecular force are similar except the 
obvious differences of the energy and length scales and 
answers the fundamental question why does the nucleus 
is approximately a collection of nucleons rather than a 
big "quark bag". Glozman and Riska argued that the 
proper effective degrees of freedom of low energy QCD 
are constituent quarks and Goldstone bosons only and 
proposed the Goldstone boson exchange model 8]. It 
also gives a good description of the baryon spectrum but 
still in its infancy for the NN interaction Q . This model 
will have the same drawback as the hybrid model in the 



understanding of the NN interaction and nuclear struc- 
ture. 

All of these models enjoy the phenomenological suc- 
cess in the description of baryon properties and even 
the baryon-baryon interactions [13] ■ The disagreement 
with the experimental data can generally be repaired by 
fine-tuning of model ingredients and parameters. How- 
ever, the more the parameters, the poorer the predic- 
tion power. In this respect QDCSM employs the least 
adjustable parameters to explain more physics related 
to the NN interaction and nuclear structure. However 
in the quantitative fit to the NN and Hyperon-Nucleon 
scattering, QDCSM is not as good as the hybrid model 
ones up to now [ll|, In order to explore if the 

QDCSM can fit the scattering data quantitatively well, 
a further extensive calculation is being done. 

This paper reports the results of the NN phase shifts 
calculated with the extended QDCSM model, where one- 
pion-exchange (OPE) with short-range cutoff is added to 
the original QDCSM. Because the effective long range 
NN interaction of the original QDCSM has been found 
to be decreased too fast due to a Gaussian wave func- 
tion used for the quark orbital motion and so the well 
established long range tt exchange is missing. This ex- 
tended QDCSM reproduces the deuteron properties very 
well 0, US]. In this calculation the tensor part of one 
gluon exchange(OGE) is also taken into account. A new 
method of calculating the phase shift is tested. Instead of 
the usual asymptotic connection to the scattering wave 
function, we first solve the equation with a zero bound- 
ary condition as one does for the bound state problem, 
by varying the boundary to check if the solution is an ex- 
ponential decreasing bound state or an oscillating scat- 
tering state, then compare the oscillating wavefunction 
with that of free particle to obtain the phase shift. This 
has been used in the lattice QCD calculation of the phase 
shift and so we check it with this otherwise can be cal- 
culated by the normal method. 
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This paper is organized as follows: in Sect. II. a brief 
introduction of QDCSM and calculation method is given, 
results and conclusions are given in Sect. HI. 
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II. HAMILTONIAN, WAVE FUNCTIONS AND 
CALCULATION METHOD 

The details about QDCSM can be found in the Refs.P, 
Here only the model Hamiltonian, wave functions 
and the necessary equations are given. The resonating- 



where all the symbols have their usual meaning. OPE 
potential and the tensor part of OGE (which is omit- 
ted before) have been added. The OPE was introduced 
to account for the long range tail of NN interaction and 
the tensor part of OGE and OPE were included to do 
the S-D channel coupling. In the OPE potential, Q(r) is 
the cutoff function, which is introduced to avoid double 
counting in the short range part. The cutoff parame- 



where fc is the channel index, Si^i = 1, • • • , n are gener- 
ating coordinates. For example, for the deuteron bound 
state calculation, we have k = 1,...,5, corresponding 
to the channels NN 5* = 1 i = 0, AA 5 = 1 L = 0, 
AA S' = 3 L = 2, NN 5 = 1 L = 2, and AA = 1 L = 
2. ij^'s are the delocalized quark orbital wave functions. 



group method (RGM), used to calculate the scattering 
phase shifts, can also be found in Refs.[l^, [l^. 

In QDCSM, the Hamiltonian for the 3-quark system is 
the same as the usual potential model and for the 6-quark 
system, it is assumed to be 



(1) 



I 

ter ro and color screening parameter /i are determined 
by the properties of deuteron. The quark-pion coupling 
constant /g^^ can be obtained from the nucleon-pion cou- 
pling constant Jnntt, which is determined by experiment. 

Combining RGM and generating coordinates formal- 
ism, the ansatz for the 6-quark system wave function can 
be written as [HI, 



(2) 



N{e) = v/l + e2 + 2ee-^'/4fc2. (3) 



62 6 
i—l i<j — l 



if i,j occur in the same baryon orbit, 
if i,j occur in different baryon orbits. 



1 + -o-^ • a J 



1 



1 1 1 

+ +3 



0{nj - tq) = 



rij < ro, 

1 otherwise, 



k i=l Lfc=0,2 ^ Q=l ,3=4 
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rj and x ^re the spin-isospin and color wave functions. 
The delocahzation parameter e is variationally deter- 
mined by the six-quark dynamics. 

With the above ansatz, the RGM equation becomes an 
algebraic equation, 



Lk 



H 



k' ,L' ,,k,Lu 



(4) 



where N'^ , H^'^^'^ are the (eq.®) wave function 
overlaps and Hamiltonian matrix elements, respectively. 
By solving the generalized eigen value problem, we ob- 
tain the eigen energies of the 6-quark systems and the 
corresponding wave functions. 

It is well known that the scattering wavefunction of 
the two-body scattering process is the spherical Bessel 
function in the asymptotic region (if there is no Coulomb 
interaction or it is neglected) but has an additional phase 
shift in comparing to the free particle ones. So the 
eigen wavefunction obtained from the eigen value prob- 
lem (eq.Q) of the scattering state should be proportional 
to the spherical Bessel function with a phase shift out- 
side the interaction range. To fix the phase shift, we 
compare the nodes of the obtained wavefunction with 
the roots of jhikR + 5l) = 0. In this way, a set of linear 
algebraic equations are obtained. The scattering energy- 
momentum and phase shift then can be obtained by solv- 
ing these equations. In practical calculation, a boundary 
condition with finite range (finite R) is always used, so 
the boundary effect should be taken into consideration 
when pick up nodes. 

Comparing to the conventional method in obtaining 
the phase shift, the new method usually takes more com- 
puting time, and doesn't work well for very low energy 
scattering because there are not enough nodes in the ob- 
tained wavefunctions. However the new mcthodhas the 
advantage that it does not need to set up an asymptotic 
boundary, where the matching conditions are used to find 
the phase shift. Especially for those cases, where to re- 
ally get an asymptotic solution is practically too difficult 
and the conventional method cannot be used. The typi- 
cal case is the lattice QCD calculation, where one always 
deals with the problem in a limited volume. This new 
method at least is the one to tackle the scattering for 
lattice QCD. It is also a convenient method to calculate 
the bound and scattering channel coupling problem. 



III. RESULTS AND CONCLUSION 

Most of the parameters in the QDCSM have been fixed 
by matching baryon properties, except for color screening 
parameter /i and the short-range cutoff tq. In this cal- 
culations, we test three values of ro- 0.6 fm, 0.8fm and 
1.0 fm. For each cutoff, fj, is determined by matching 
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FIG. 1: N — N phase shifts versus laboratory kinetic energy 
for the channel of ^''^Si. 



the mass of the deuteron. To save space, here we only 
report the results for the best one set (ro — 0.8fm) [1J| 
For comparison, the results without OPE are also given. 
The model parameters aregiven in Table I. 

Table I. Model parameters 





without OPE 


with OPE 


ro (fm) 




0.8 


m (MeV) 


313 


313 


6(fm) 


0.6034 


0.6015 




1.543 


1.558 


ac (MeV fm-2) 


25.132 


25.135 


fi (fm~^) 


1.0 


0.85 


(MeV) 




138.04 



The following NN phase shifts, ^.i^^, 
Di, ^'^D^, have been calculated. 



The symbol 



1,3 

2/+i,25+i^ ^ is used in the paper and all the experimental 
data of NN phase shifts are taken from Ref . [ij, Ull . For 
the case of ^^'^S'l, ^'^Di, coupling between S- and D-wave 
is considered. 

For ^''^S'l (FIG.l), the phase shifts in the single chan- 
nel calculation decrease when OPE is added due to the 
smaller fi used. OPE provided additional attraction and 
stronger S-D coupling, the color screening fj, has to be re- 
duced for compensation to reproduce the deuteron bind- 
ing energy. Taking into account the effect of S-D mixing, 
the phase shifts rise and the results approach the ex- 
periment data again. The remaining discrepancy might 
be resolved by additional channel coupling and this has 
happened in the bound sate calculation of deuteron. 

The ^'^Sq channel has a similar results which is not 
given to save the space. 

For the case of ^'^Di (FIG. 2), OPE improves the phase 
shifts, especially when S-D mixing is considered. The 
S-D mixing pushes the phase shifts down, so does the 
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FIG. 2: Same as FIG.l for the channel of '•''L'l. 
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FIG. 3: Same as FIG.l for the channel of ^"^Da- 



decreasing of color screening parameter due to the ad- 
ditional attraction contributed by OPE. The total effect 
makes the results match the experimental data very well. 

As for ^'^D^ (FIG. 3), the results match the experi- 
mental data as well. As in the ^''^Di case, the calculated 
phase shifts go down to the experimental data due to the 
decreasing ofthe color screening parameter after adding 



the OPE. 

These results show that QDCSM is possible to de- 
scribe the NN scattering data quantitatively by includ- 
ing the long range OPE and channel coupling. Tak- 
ing into account the following two results together, that 
the deuteron properties have been reproduced very well 
by the same extended QDCSM 0, [H and the effective 
baryon-baryon interactions of the three constituent quark 
models have been shown to be similar . it is fair to 
conclude that the short and intermediate range NN in- 
teractions, which were attributed to the a (or two tt) 
and heavier meson (p, w, etc.) exchanges in the Gold- 
stone boson exchange and hybrid model, are successfully 
described by the quark delocalization and color screening 
effects. 

Quark delocalization is exactly the same as the elec- 
tron delocalization which is responsible for the molecular 
bond. The nuclear force and molecular force share the 
same delocalization mechanism and therefore (at least 
in our opinion) appear to be similar. As for the color 
screening, it is a phenomenology to describe what has 
been missing in the two body confinement. Some non- 
linear QCD interactions, such as the three gluon interac- 
tion, three body instanton interaction j_lS||, certainly can 
not be included in the two body confinement. However a 
real QCD verification of the model Hamiltonian assumed 
in QDCSM is still an expectation. 

For the quantitative fit to the scattering data, there 
are still disagreement remains to be resolved in QDCSM. 
The spin-obit interaction of the effective one gluon ex- 
change should be added to understand the P-wave and 
higher partial wave phase shifts. The A^A, NT, scatter- 
ing should be reanalysed with the Extended QDCSM. 
The NE scattering should be predicted before the ex- 
perimental measurement and it will be a good check of 
different constituent quark models. 

Our experience on the new method of calculating the 
phase shifts is preliminary. Further quantitative check 
with the usual method within larger energy region and 
more channel coupling is needed. Our message to the lat- 
tice QCD phase shift calculation is, it is a unified method 
to deal with the bound state and scattering but should 
be used with caution especially in the very low energy 
region and quantitative reliable results being expected. 
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